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SUMMARY 

The  present  work  is  a  further  step  towards  optimization  of  stiffened  cylindrical  shells*  Rings 
and  stringers  of  non-uniform  cross-section  are  analysed  as  a  means  of  obtaining  higher  structural  effi¬ 
ciencies.  For  lateral  pressure  loading^ings  with  non-uniform  height  or  width  are  compared.  The  gain 
in  general  instability  pressure  relative  to  the  uniformly  stiffened  she’’,  of  the  same  weight  amounts  to 
70  —  80%.  For  prescribed  loading.  10  —  20%  weight  savings  are  obtained.  For  axially  compressed  cyl¬ 
inders  stringers  with  sinusoidal  and  linear  height  variation  are  studied  and  compared.  Gains  in  load 
and  weight  savings  of  30%  and  10%,  respectively,  are  obtained,  and  the  sinusoidal  variation  is  found 
to  be  superior  to  the  linear  variation. 

Finally,  cylinders  under  hydrostatic  pressure  stiffened  by  a  combination  or  uniform  stringers 
and  non-uniform  rings  are  analysed.  Different  configurations  are  checked  in  order  to  obtain  tbe  highest 
possible  efficiencies. 
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1.  INTRODUCTION 

Experimental. and  theoretical  iaveatigationa  of  uniformly  stiffened  cylindrical  shells  with  close¬ 
ly  spaced  heavy  stiffeners,  [1]  to  [6]f  have  shown  that  for  axial  compression!  lateral  pressure  or  torsion, 
the  deflection  shape  during  general  instability  consists  of  one  half  wave  in  the  axial  direction.  Excep¬ 
tions  to  this  rule  are  the  cases  of  ring  stiffened  shells  under  axial  compression  or  hydrostatic  pressure, 
in  which  buckling  occurs  with  many  waves,  see  [1]  and  [3], 

One  approach  to  higher  structural  efficieacies  is  to  use  non-uniform  stiffeners,  the  intensity  of 
which  has  a  maximum  at  the  mid-length  of  the  shell,  •  A  similar  approach  used  earlier  in  optimum  design 
of  columns  [7]  and  stiffened  plates  [8]  achieved  considerable  weight  savings  by  variation  of  bending 
stiffness  along  the  length  or  width  of  the  structure. 

In  shells,  this  variation  in  stiffness  can  be  obtained  in  two  ways :  with  uniform  stiffeners  and 
varying  spacing,  or  with  uniform  spacing  and  varying  stiffener  cross-section.  Obviously,  a  combination 
of  both  is  also  possible.  One  may  note  that  on  account  of  edge  effects  both  unstiffened  and  uniformly 
stiffened  cylindrical  shells  are  actually  structures  of  non-uniform  stiffness  in  the  axial  direction.  Non- 
uniform  distribution  of  stiffener  material  can  compensate  for  the  edge  effects  and  stiffen  the  shell  at  the 
weakest  portion.  Hence  U  represents  on  approach  to  optimization. 

For  a  cylindrical  shell  unequal  spacing  yields  *aub-#hella"  (between  rings)  of  non-uniform 
buckling  strength.  The  sub-shell  next  to  the  bulkheads,  where  the  largest  spacing  occurs  is  the  most 
prone  to  buckling.  The  increase  ia  general  instability  will  therefore  be  overshadowed  by  s  decrease  in 
local  stability.  This  argument  has  to  be  born  in  mind  all  along  the  analysis  aince  optimization  is  here 
mainly  concerned  with  general  instability.  In  contrast  to  cylindrical  shells,  conical  shells  are  examples 
of  structures  where  non-uniform  ring  spacing  ia  very  efficient.  In  [9]  and  [101  an  optimum  weight  analy- 
aia  yielded  weight  savings  of  more  than  25%  for  non-uniform  ring  spacing. 

Variation  of  stiffener  cross-section  is  the  alternative  method  of  stiffness  variation  along  the 
length.  In  cylindrical  shells,  such  an  approach  to  general  instability  optimization  does  not  influence 
the  instability  of  the  sub-shells.  It  is  assumed  that  the  influence  of  stiffener  rigidity  on  local  shell  in¬ 
stability  is  small.  Since  local  stiffener  instability  ia  rarely  a  design  criterion,  the  change  in  stiffener 
buckling  strength  is  not  considered. 


The  theoretical  analysis  was  performed  on  st  gera  and  rings  of  rectangular  cross  section.  The 
height  or  width  of  the  stiffeners  may  be  changed  aloug  the  generator.  In  some  cases  a  combination  of 
both  variations  yields  the  most  efficient  structure.  '  Some  typical  variations  of  stiffeners  are  shown  in 
Fig.  2.  Note  that  in  practical  designs  only  part  of  the  width  or  height  of  the  stiffeners  is  varied.  Height 
variation  is  clearly  more  effective  since  the  moment  of  inertia  of  the  stiffener,  which  determines  its 
main  contribution  to  the  stiffness  of  the  shell,  varies  according  to  the  third  power  of  the  height.  The 
maximum  permissible  height  is,  however,  fixed  by  practical  considerations,  such  as  for  example: 

1.  The  thickness  of  commercially  available  plates  for  stiffened  cylinders  of  integral  construction* 

2.  The  maximum  feasable  slenderness  ratios. 

These  constraints  are  sometimes  so  severe  that  the  relatively  inefficient  width-variation  is  pref¬ 
erable.  It  may  be  mentioned  that  the  machining  difficulties  which  arise  in  the  case  of  height  variation 
do  not  appear  in  the  case  of  width  variation. 

Both  types  of  cross-section  variation  and  different  laws  for  this  variation  (linear,  sinusoidal, 
see  Fig.  2)  are  considered.  The  type  of  stiffeners  to  be  used  (rings,  stringers  or  combination  of  both) 
depends  mainly  on  the  external  loading.  For  lateral  pressure  and  torsion,  rings  alone  are  the  most  ef¬ 
ficient,  see  [1]  and  [2]  while  for  axial  compression  a  combination  of  rings  and  stringers  is  the  best, 

[10]  and  [18]. 


2.  THEORETICAL  ANALYSIS 

2.1,  DIFFERENTIAL  EQUATIONS  AND  BOUNDARY  CONDITIONS. 

The  analyeis  starts  with  the  same  main  assumptions  as  [11].  For  the  sake  of  clarity  they  are 
repeated  here. 

a)  The  stiffeners  are  distributed,  or  "smeared",  over  the  whole  surface  of  the  shell.  A  more  rigor- 
ous  analysis  that  takes  the  discretness  of  stiffeners  into  account  is  given  in  [12]  and  [13].  it  is  shown 
there  that  unless  the  number  of  rings  is  less  than  3-4  the  difference  between  the  buckling  load  predicted 
by  'distributed"  and  "discrete"  theories  is  negligeabie. 

b)  The  normal  strains  tx  and  vary  linearly  in  the  stiffener  as  well  as  in  the  sheet.  The  nor¬ 
mal  strains  in  the  stiffener  and  in  the  sheet  are  equal  at  their  point  of  contact. 
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c)  The  stiffeners  do  not  transmit  shear*  The  shear  membrane  force  is  carried  entirely  by  the 

sheet. 

d)  The  torsional  rigidity  of  the  stiffener  cross-section  is  added  to  that  of  the-sheet  (the  actual  in¬ 
crease  in  torsional  rigidity  is  larger  that  that  assaued). 

e)  The  strain  displacement  relations  are  similar  to  those  used  by  Donnell  in  [14]. 

The  equations  of  equilibrium  for  buckling  under  axisymmetrical  loading  are 


N*.« +  “  0 


Nij+Nj  » 0 


M«.«  +  +  +  RNtf  +  RN*0  w.*«  +  RN^0  w.^  *  0 

where  z,  x  and  <*>  are  non-dimensional  co-ordinates  (see  Fig.  1). 

The  stress-strain  relations  differ  for  shell  and  stiffeners  and  are 


(1) 


o,(z)  =  (E/(l-„2)Hfx  +  Vt^  -  z  (kx  +  VK^i\  ^ 


ff^(*)  "  [E/(l-i'2)][f0  +  r  zGc^  + 
(z)  -  E  (<_  -  z  *  ) 


V*)  *  E  ~z  V 


>  -  in  the  shell 


L  -  in  the  stiffeners 


(2) 


Note  that  in  Eqa.  (2)  the  middle  surface  of  shell  is  used  as  reference  line.  Hence,  the  for ce3 
and  moments  acting  on  an  element  become 


NX  =  (Eh/(W2)]i;*x(l  +p*(x)]  +  ^-Kxx*  U)i 
N^  =  (Eh/(1-^2)]|(^[1  +  P*(x)]  +  V(t  -  1 

[Eb/2  (1  +  »/)]  yx^ 

M,  =  -  (D/R)Ui[l  +  qjj  (x)  ]  +  ~  ixC[  (*)  1 

%  -  -(1T/B;U^[1  +  t/*2  (x)  ]  T  V«K  Tf'fig  (x)l 


(3) 


M«<£  =  D/R  [1  -  V  +  (x)  ] 

M^--(D/R)  [I-**,*  («)]«„* 


(3) 


The  curvatures  in  Eqs.  (2)  and  (3)  are  non-dimensional,  the  physical  parameters  having  been 
divided  by  R.  p|(x) ,  fx*(x) ,  (x) ,  :;*2  (x) ,  »/*j  (x)  and  (x)  are  the  changes  in  stiffness  due 

to  non-uniform  stringers  and  rings  while  y*{x)>  X  J(x),  £*(x)  and  C Jd)  represent  similar  contribu¬ 
tions  due  to  the  eccentricity  of  stiffeners  (positive  for  inside  stiffeners  and  negative  for  outside  ones). 
Their  definitions  in  terms  of  stiffener  geometry  for  any  type  of  cross-section  and  for  the  particular  case 
of  rectangular  cross-section  are  given  in  Appendix  A.  Note  that  for  constant  stiffener  cross-section 
Eqs.  (3)  degenerate  to  Eqs.  (5)  and  (6)  of  [111. 

The  strains  and  curvatures  are  defined,  as  in  [  14],  by: 


V  "  V.*~W 

vv., 

K_  «-  W 


(4) 


.«« 


K<t>  ’ 


W.x* 

Substitution  of  Eqs.  (3)  and  (4)  into  Eqs.  (1)  yields  the  equations  of  equilibrium  in  terms  of 
displacements 

*  l(1  -  ^>/21  +  K1  +  v)/2]  v.,<£  +  K*'*? (x)1.«  ‘  [w.«  (x)1.«  "  0 

\<H>  +  td-M  V  „  +  [(l  +  v)/2]  u  +  (v  ^  -  w>0)  pj(x)  -  X*  (x)  =  0 

y4  W  +  ’01  (x)l,„  “  K*  ^2  to  +  2  -  V.W  ^  (x>  * 

+  +  »??..  (*)ll+12(R/h)  Kw-v  j)11+P?(x)1-«/u  J+(pR  /D)w  +(PR/2nD)w  -0  (5) 

»xyr  •  *  *•  »*  »—  tyr  #xs 


.t 


}  .  . 
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The  prcbuckling  membrane  stresses  for  a  combination  of  axial  and  lateral  loadings  are : 

N  •>  -  (P /2  nr  R) 

(e 

^0  =  -pR 

Equations  (5)  are  a  set  of  three  linear  homogenous  partial  differential  equations  with  varying 
coefficients.  The  usual  simply  supported  boundary  conditions  are  assumed: 


Mx  =  N,  =  0 


w  =  v  =0 


x  =  0  and  L/R 


Two  methods  of  solution  are  proposed.  In  the  first,  a  straightforward  Calerkin  procedure  is  ap¬ 
plied  to  the  three  equations  transforming  them  ic*,o  infinite  sets  of  linear  algebraic  equations.  The 
second  employs  an  approximate  procedure  which  "corrects"  the  first  two  stability  equations  and  then 
solves  the  third  one  by  a  Calerkin  procedure.  Numerical  comparisons  between  both  methods  are  carried 


2.2.  SOLUTION  BY  GALERKIN'S  METHOD 

The  variational  equivalent  of  Eqs.  (5)  is  written  as 
L/R  2  rr 

(  f  lu#„  +  Kl-pVSlu^  ,£  +  [(l  +  p)/2]  v  -  j/w^  +  [u(J(fiJ(x)]>jt  -  [w^xJMljSu  dx  d«/>  =  0 

L/R  2 17 

f  i  |v,<W>+((1~l/)/21v.«  +  £(1+l/)/2K*<i  “  w,<£  +  (\<M>  -  w,<t)  **2  (x)  -  dxd«B  0 

0  0 
L/R  277 

f  W+fW.xx'lo*l(x)U~[u.*^t(x)U"W,<Mi'702(x)+2(w.^-V,<^)^2(,t>  + 

0  0 

+  +  7*2W1I  +  12  (R/h)  |(w  -  V>0)  [1  +  P*(x)i  -  vu  J  +  (pR  /D)w^  + 

+  (PR/2  77  D)  w  a  ^  8w  dx  di/i  •»  0  (8) 

The  displacement  components  are  expanded  into  Fourier  series  in  the  axial  direction 
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u  *  <»n  t  <j>  2  A  cos  n  fl  x 

n*l  ° 

V  -  costal  Boflin  n£x 

w  -  «»n  t  <6  2  C  *b  n  »x 

Sene»  (9)  fulfil  the  boundary  conditions,  Eqs.  (7). 

A  permissible  variation  of  Eqs.  (9)  is 

=  8*n  t  <6  cos  m0  x  S  A 

IB 

Svm  *  cos  t  <f>  sin  mfi  x  8  B 

m 

Swm  =  «'»  *  *  »i«  m/3  x  5  C 

IB 

Substitution  of  Eqs.  (9)  and  (10)  into  Eqs.  (8)  yield. 

•O 

Fa  (».*)+  Bn  Fe  (n,m)  +  CnFc  (a,B)i  sA,  .  o 

oo 

«Ii  tA0GA(nfm)  +  BnGB(„,m)  +  CBGc(n,«)]  SB,  -0 

«l  £AA  <*.»>  +  BbHb  (n.m)  +  CaHc  (n.m)I  SC,  -  0 
where  tbe  F'.,  G's  and  H’s  are  defined  a. 

Fa(«.»>  -  lo8^2  +  [(1  -  v)/2)  t2  I  8am  +  amB*lj(n,m) 
FB(n,a)  *  CA{nim)  *  ((1  -  j/}/2]  tn;3  8 

Fc(“’m)  "  Ha (m’n>  «  ^BSBm-a2»BJl8(..») 

C3  (n.m)  -  It2  +  [{1  ~  j/)/2]  n2  /92|  S  +{2K  f-  ~) 
Gc(n,m)  «  Hb  (u,m)  .  1  [Sam  +  Kj  (n,m))  -  t3  (n>la) 


(10) 


(11) 


(12) 


(13) 
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Hc(n,ra)  =  |1  +  [t2  +  n2^2)2/12  (R/b)2]-[(A  t2  +  An2  $2/2)/12(R/h)2]  I  8  -2t2K,(n,m)  + 

v  n  m  z 

+  Kj  (n.ra)  +  t4K3(n,m)  +  i  fn2 m2  ft*  In  (n,m)  +  t2nm  £2  ^  (n,m)]/)2(R/b)2|  (14) 


Ap  and  A  are  non-dimensional  lateral  and  axial  loading  parameters  respectively.  I.(n.ra)  and  K.  (n.m) 
are  weighted  integrals  of  the  functions  describing  the  variation  of  stiffener  cross-section.  They  are 
defined  by 

(rr/fi) 

K  (n,m)  =  (2 il/rr)  f  f  (x)  sin  n  Bx  sin  m/Sxdx  \~  1,2,3 

1  o  ' 

f,(x)  =  ,r*(x)  f2(x)  =  v|(x)  f3(x)  -  » )*2(x) 

in/fi) 

I,  mi. in)  -  (2/m ;r)  f  t/r*(x)  sin  n  /3x)  g  cos  m  f3x  dx 
(tr// 3) 

I2(n,rn)  -  (2/mir)  £  J  (x)  sin  n/?x]  xcos  m/3xdx 
(ir/f  J) 

1,  (n.m)  =-(2/m 2nft)  f  (^j(x)  sin  n  /9x]  a  sin  m  /9x  dx 
(rr//3) 

l^U.m)  =-{2/m n)  f  (x)  +  i]f2{x))  cos  n  /3  x]  x  sinm/3xdx 


In  all  derivations  it  was  tacitly  assumed  that  the  functional  variation  of  the  stiffener  cross  sec¬ 
tion  is  a  function  of  Class  2  (continuous  up  to  derivatives  of  the  second  order)  otherwise  some  of  the 
integrals  of  Eqs.  (15)  are  meaningless.  Note,  that  the  linear  variation  (see  Fig.  2)  exhibits  a  disconti¬ 
nuity  in  slope  at  x  =  L/2R.  This  case  and  ones  with  a  similar  type  of  discontinuity,  must  be  handled 
by  a  different  method.  Eq3.  (5)  which  include  derivatives  of  the  functional  variation,  have  different 
meanings  when  x  is  on  the  right  or  when  it  is  on  the  left  of  the  discontinuity,  A  detailed  analysis  for 
such  cases,  that  assumes  symmetry  of  stiffener  distribution  relative  to  x  =  L/2R,  is  given  in  Appendix 
B. 


After  repeated  integration  by  parts  of  Eqs.  (15),  the  final  results  are 


K .  (n.m) 
1  (n.m) 

K 


(2fi/rr)[l  +  (-  l)0+<n] 


(sr/2/S) 

f  f.(x)  sin  n  fix  sin  m  j9x  dx 
o  * 


,  4  r.  *-m  . 

{2fi/n)[  1  ¥  (-  1)  )  (  g.  (x)  sin  n  /3x  sin  m  /Sx  dx 

o  J 


1,  2,  3 


(16) 


I 
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a+m  (ff/2/3) 

I4(n,ra)  *  (2/5/»r)[l  +  (-1)  ]  f  U?*i  (x)  +  »;*2  (x)]  cos  n /dx  cos  m /9x  dx 

gj(x)  =  pj(x)  g2(x)  =  x*(x)  g3(x)  =  j7q  3  (x)  (16) 

Eqs.  (16)  apply  to  variation  of  stiffener  cross  section  with  continuous  or  discontinuous  slopes. 

If  one  now  returns  to  Eqs.  (11)  and  uses  the  arbitrariness  of  the  variations  8k  ,  SB  ,  SC  one 

m  m  m 

obtains  three  sets  of  linear  algebraic  equations  which  can  be  written  in  compact  form  as 


Consideration  of  only  the  first  N  terms  of  series  (9)  yields  stability  determinants  Eqs.  (17)  of 
order  3N .  Each  element  of  Eq.  (17)  is  an  N  XN  square  matrix  of  which  only  l  Hc  1  contains  the  ex¬ 
ternal  load  along  its  diagonal.  The  lowest  eigenvalue  of  Eq.  (17)  yields  the  critical  load  for  general 
instability.  The  integral  value  of  t  (the  number  of  circumferential  waves)  which  minimizes  the  critical 
load  must  be  used  in  compulations. 

According  to  Eqs.  (12).  (13),  and  (14) 

fGAJ  -  fFBJ  [Rb]  -  [Gc]  [Ha ]  -  [Fc]T  (18) 

T 

In  Eqs.  (18).  [  ]  represents  the  transpose  of  the  mntrix  and  [Hg]  and  [Gc]  are  symmetric. 

The  stability  matrix.  Eq.  (17),  is  therefore  symmetric  as  would  be  expected  from  Maxwell's  re¬ 
ciprocal  theorem.  This  symmetry  is  of  great  help  during  the  numerical  computations.  Integrals  Lfn.m) 
and  Kj(n.m),  see  Eqs.  (16'  e  identically  zero  for  odd  values  of  the  integer  (n  +  m).  Under  these 
conditions  the  stability  detek.  tant  of  Eq.  (17)  can  be  resolved  into  two  sub-determinants,  one  of  the 
even  components  and  one  of  the  odd  components  of  the  displacements.  The  two  sub-determinants  re¬ 
present  symmetric  and  antisymmetric  buckling  modes.  The  critical  load  has  to  be  computed  for  both  pat¬ 
terns,  and  the  one  yielding  the  lower  load  is  the  one  to  be  considered. 

As  mentioned  earlier,  the  unknown  in  Eq.  (17)  appears  only  in  lHr  ).  This  feature  is  used  to 
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develop  a  systematic  solution  of  Eq.  (17)  whose  outline  is: 

1*  Application  of  Gauss’s  algorithm  to  reduce  the  order  of  matrix  from  3N  to  N 

2.  Use  of  an  iterative  procedure  for  calculation  of  the  lowest  eigenvalue  and  eigenvector  of  the 

reduced  matrix. 

Details  are  given  in  Appendix  C. 


2.3.  THE  METHOD  OF  “CORRECTING  COEFFICIENTS** 


For  a  uniformly  stiffened  cylindrical  shell  the  differential  equations  of  equilibrium,  Eqs.  (5),  in 
the  presence  of  the  boundary  conditions  for  classical  simple  supports  Eq.  (7>  have  a  closed  form  solu¬ 
tion,  see  [3]  or  [111.  A  non-uniform  distributite  of  stiffeners  increases  the  stiffeness  of  the  middle 
part  of  the  shell  relative  to  that  at  the  edges.  Daring  buckling,  the  deflection  curve,  which  in  the  case 
of  constant  stiffeners  is  a  half  sine  wave  in  the  axial  direction,  will  be  of  similar  shape  with  a  pronounced 
flatness  at  the  middle  of  the  shell.  If  one  considers  the  expressions  of  displacements,  Eqs.  (9),  this 
means  that  the  first  term  of  the  series  is  dominant,  and  the  additional  terms  correct  the  basic  mode  only 
slightly  without  altering  its  general  features.  This  is  the  main  assumption  of  the  method  of  "correcting 
coefficients"  used  earlier  [15)  in  connection  with  general  instability  of  stiffened  conical  shells.  In  the 
application  of  the  method,  the  varying  coefficients  of  the  first  two  differential  equations,  Eqs.  (5),  are 
changed  to  constant  ones.  Substitution  of  series  (9)  into  the  resulting  two  "corrected"  equations  gener¬ 
ates  an  infinite  set  of  uncoupled  linear  algebraic  equations.  Hence  only  the  general  harmonic,  n,  has 
to  be  considered.  The  third  stability  equation  with  increments  from  the  first  and  the  second  is  solved 
by  a  Galerkin  procedure. 

The  variational  equivalent  of  differential  equations  (1)  is 


2rr  (L/R) 


»v  .  ;  j  l|N,.,  .  N.^1  4.  ♦  [N^  .  S,  ♦  [M,.„  .  -MtsW 

-  *  RN*  *  RN.«  ”...  *  (WI01  *  N.*« 


V  - 


-  (M,/BXSwifi  +  [(Ml#lt/R)  -  +  N,o  J 5w!  R  ^  -  0 


(19) 


This  equivalent  is  the  variation  of  the  generalist  .!  potential  after  two  integrations  by  parts. 
Since  the  variations  are  arbitrary,  the  differential  equations  (1)  and  boundary  conditions  (7)  can  ready- 
ly  be  obtained  from  Eqs.  (19). 

To  the  membrane  forces,  defined  by  Eqs.  (3),  a  set  of  terms,  whose  sum  is  2ero  is  added.  Then 


Nx  -  [Eh/(1  -  »<2)1  lo^(l  +  kj)  »  t/(v^-wr)-w>](!tk2l  +  [Eh/(l-^2)]lutX[gj-k1l- 

~w. k2]{  “  +  N„  (20) 

where  kj  and  k2  are  constants.  The  second  term  of  the  last  expression  in  Eq.  (20),  Nx#,  is  called 
"Error  Load** 'Similarly  the  membrane  force  is  transformed  to 

*  [Eh/(!  —i'2)]  Uv^  -  w)  (1  +  kj)  +  vu,x  “  ^4  + 

+  [Eh/(1  -  i/2)]  l(v^-  w)  [n*{x)  -  k3l  -  w  _  kt]  }  (21) 

Substitution  of  Eqa.  (20)  and  (21)  into  Eq.  (19)  yields 

2»r  (yaj 

l  [  \lK*  +  N^.^3  5“  +  +  N*0,«3  5V  *  jN*.,x  Su  +  N0..05v  + 

+  +  +  ~  +  RN0  +  RN0O  ",<£<£  +  RN«o  W,.J 

(Sw/R)}}  dxd<£  -  0  (22) 


since  the  line  integral  vanishes  on  account  of  boundary  conditions,  Eqs.  (7),  which  are  satisfied  by  the 
assumed  displacement  aeries. 

The  corrected  first  two  stability  equations  are  then 


N 


N 


N 


6 


»* 


0 


(23) 


Assumption  of  a  solution  of  the  type 
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u  =  Ab  sin  t  $  cos  n$x 

v  *  Bo  cos  t  <£  sin  u(3x  (24) 

w  -  C(  sio  t  <f>  sin  n/3x 

and  use  of  Eqs.  (20)  and  (21)  yields  two  homogenous  algebraic  equations  with  three  unknowns  An,  Bn 
and  Cn,  which  are  related  to  each  other  through 

A  =  a  c  a=(D,/Dn) 

a  a  a  a  '  Is'  Ob 

(25) 

E.  ■  ‘.S  l.  - 

where 

D0n  -  ((l-(U2)(l+k  (l+k,)(l +k,)-vlnl^tJ +[(l-v)/2](l+ki)n4/34 

D i „  =  Ml  +  v)/2]  ni9t4k4  +  (1  +  k3)  [k,  ns  f?  +  [(l-v)/2]  a/3l  t2  +  [(l-„)/2)  n3  [k2  n2  $  -  v] 

DJn  -  {(1  -  v)/2\  k4 15  ♦  KUkj)n2  /3*k4  -  ((l-v)/2](l+ks)  It3  +  $\  [Al  +  v)/2]  - 

(1  +  k.)(l  +  k,)  ln2/32  -[(1  +  vV2]  k2n4/34)  t  (25a) 


The  third  equation,  with  the  "Error  Loads", 

2  rr  L./R 


f  f  IN  5u+N,  ,5v+(M 

o  o  **•*  9»>«P  *'* 


M 


+  RNqSo  W,<^  +  RN.0  ".a*1  (Sw/  R)j  ix  ** 


(26) 


i«  solved  by  a  Galerkin  procedure.  Substitution  of  Eqs.  (25)  into  tho  expressions  of  the  displacements, 
Eqs.  (9),  and  their  variations,  Eqs.  (10),  yields 


u 

v 


sin  t  -fr 
cos  t  <p 


OC 

2 


a  C  coa  sBx 

b  n 

bnCo  sin  n^x 


5um  «=  a^  sin  t$  cos  m  /3x  8Cm 

8v  =>  b  cos  t<£  sin  m  fix  5C  (27) 

IB  W  'in 
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w  -  lb  t f  C  ain  a/3x  5wm  -  .in  t  <6  sin  m  0  x  SC 

Substitution  of  Eqa.  (20),  (21)  and  (3)  in  Eq.  (26)  yields 

2ff(L/It) 

0  “  f  f  #Kl5u  +  (Eh/(l-v2)l 

* (v,^  *  ",*>  I<*2  (x>  ~  k3]  ~ Uf  2*  (*>  -  k4ll  Sv  -(D/R2)  |74 w  +  (w  „  ^  (x)  ]  „ 

-Kx  <r<x>  u  ♦*.**«:* (x> + »■.*- w  «  w  ♦  K^rt  <x> + 

+  *7*2 ^X^,»  +  12(R/b)  [(v^-  w )L\+n*(x))  -  vu>xJ  +  kp  w,<fxf>  +  A(w  x,/2)!5w}dxd<^ 
Use  of  Eqs.  (27)  and  the  definitions  of  the  weighted  integrals,  Eqs.  (16),  yields 
S  T(n,  m)  C  =  0  m  =  1,  2 . . 

»«i 

where 

T(q,ni)  ^  (Q(n,m)  +  1  [(n202  +  t2)/ 12(R/h)2]  +  1  +  b# t  +  »a/3«B -[(An2/j2/2)  + 

+  Att2]/12(R/h)2 }  Smm  +  l[x2  m2  0*  I,  (n,  m)  +  n  m  0  t2 I4  (n,  m)  +  t«  K,  (n,  m)]/12(R/h)2|- 
—  n  m2  0*  a%  Ij  (n,  n)  -  t2(2  +  b>t)KJ(n,m)  -r  (l+b^Kj  (n,m)) 

Q(n.  m)  the  work  done  by  the  error  loads,  is  defined  by 

Q(n,m)  =  I*. (Ij  (n, m)-kj  5>-]-n0[l2(n,nO-k25nB]}  a0  m/3  a^  + 

♦  Ml  +  b,  t.)  [K,  (a,  at)  -  k,  6,  J  -t2  [K ,  (s,  at)  -  k4  8,  J 1 1  b. 


(27) 


(28) 


(29) 


(30) 


A  finite  number  of  terms,  N,  in  Eq.  (29)  yields  a  set  of  N  linear  homogenous  equations.  For 
reasons  similar  to  those  presented  in  the  previous  section,  the  determinant  can  be  resolved  into  nub- 
detenaiuants  of  the  even  and  odd  components  respectively.  Again  the  value  of  t  that  minimizes  the 
critical  load  has  to  be  used. 
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2.4.  EVALUATION  OF  THE  CORRECTING  COEFFICIENTS 


The  accuracy  of  the  method  in  predicting  general  instability  largely  depends  on  the  manner  in 
which  the  correcting  coefficients  are  computed*  In  the  present  solution  they  are  evaluated  by  equating 
to  zero  the  work  of  the  error  loads,  and  ,  of  the  first  two  stability  equations.  This  work  is 
computed  for  different  permissible  displacements  in  the  same  modus  operand!  applied  in  Galerkin's 
method.  It  should  be  mentioned  however  that  due  to  the  limited  number  of  unknowna  (k, ,  k2,k3  and 
k4)  the  virtual  work  can  be  made  to  vanish  for  a  maximum  of  four  different  displacements. 

As  a  first  approximation  the  first  terms,  n  =  1,  of  series  (27)  are  used  as  permissible  displace' 
meats.  Bearing  in  mind  the  discussion  at  the  beginning  of  section  2.3,  this  choice  is  very  promising, 
since  n  =  1  is  the  basis  mode  of  the  deflection  shape  during  general  instability. 

The  error  work  expression,  Eq.  (30),  then  yields 

k,  =  1,(1,  1)  k3  =K4(1.1) 

(31) 

k,  =  I2(l,l>  k4  =K2(1,  1) 

i  here  are  cases,  however,  where  the  first  term  of  series  (27)  is  not  dominant,  e.g.  the  case  of 
c  .i.  ;5  stiffened  shell  underaxial  compression,  then  Eqs.  (31)  will  yield  very  unconservative  buckling 

’'uUS. 


A  mor"  systematic,  though  less  simple  way  of  computing  the  correcting  coefficients  is  to  use 
-he  actual  displacements  at  buckling  in  the  computation  of  the  virtual  work.  'In  other  words,  not  the 
first  terms  of  Eqs.  (27)  but  the  whole  series  are  used  as  permissible  displacements. 

1  hen 


N 

i  n  m  a  a  C  C  1,  (n,m) 

>  D  m  D  Ml  1 

n  ,  m  “  i 


N 

V 

n  =  1 


2  2 
n  a. 


c 


2 

o 


to  a,  CJtl,]  In  an  C„l 
In  au  Cnj  In  afl  C J 


k., 

i. 


-  n2  m  C  C  l,(n,m)  a 

•  u  m  i  m 

n ,  m  ^  1 


N 
2 
n  “  1 


i 


a 

n 


la2CoHl2]  In  ao  CB! 


in2  CnJ  In  an  Cnl 


(32) 
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N 

2  (l  +  bt)C  K.(n,m)b  C  Kl  +  b  „t)yJ  iCK.llb  Cl 

_  _ _ \  on*  ram  q  q  inn 

J,U  +  b.l>C>„  lU  +  b.OC.lll.C.I 

£  C,  „)b„C.  IC.llK.Ilb.  C.l 

I  CJb  iC  1  |b  C  i  (32) 

n  =  l 

whore  1  | ,  1  |  are  shorthand  notation  for  row  and  column  matrices. 

It  should  be  pointed  out,  however,  that  preliminary  knowledge  of  the  eigenvector  is  needed  in 
Eqs.  (32).  Hence  an  iteration  procedure  is  required.  As  a  first  approximation  one  uses  Eqs.  (31).  The 
computed  eigen-vector  is  then  used  in  Eqs.  (32).  The  process  is  repeated  until  reasonable  convergence 
of  the  solution  is  achieved. 

The  purpose  of  the  method  of  correcting  coefficients  was  to  provide  an  easier,  though  less  ac¬ 
curate,  method  for  solving  the  stability  equations  (5).  It  lead  to  an  iterative  process.  In  comparing 
computations  for  both  methods  it  is  obvious  tbnt,  in  the  present  case,  the  loss  in  accuracy  in  the  sec¬ 
ond  method  outweighs  the  numerical  difficulties  of  the  first.  However,  this  is  not  always  so,  since  for 
other  types  of  shells  (conical)  or  loads  (torsion)  a  straightforward  Calerkin  or  Rayleigh-Ritz  procedure 
might  be  very  cumbersome. 

In  the  present  work  both  methods  were  used  for  solving  a  typical  problem.  Comparison  of  the 
results  yielded  important  information  as  to  the  accuracy  of  the  approximate  method  of  "correcting  coef¬ 
ficients  ”  and  its  computional  t;r„e  saving. 

3.  NUMERICAL  RESULTS  AND  DISCUSSION 

3,1.  INTRODUCTION 

The  theoretical  analysis  considered  a  general  non-uniform  distribution  of  stiffeners,  symmetric¬ 
al  with  respect  to  the  mid-length  of  the  shell.  In  order  to  obtain  nn  estimate  of  tbe  possible  gains  in 
load  for  a  constant  structure  weight  or  of  the  weight  saving  for  a  constant  load,  one  roust  now  solve 
the  problem  numerically  for  typical  laws  of  stiffness  variation.  The  cases  considered  here  are : 
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1)  Cylindrical  shells  subjected  to  lateral  pressure  and  stiffened  by  rings  of  non-uniform  width  cr 
height. 

2)  Axially  compressed  cylindrical  shells  stiffened  by  stringers  with  sinusoidal  or  linear  height 
variation. 

3)  Cylindrical  shells  subjected  to  hydrostatic  pressure  and  stiffened  by  a  combination  of  uniform 
stringers  and  rings  of  non-uniform  height. 

All  cases  were  solved  by  the  more  rigorous  Galerkin  procedure.  Case  2  was  also  solved  by 
the  approximate  method  of  "correcting  coefficients".  The  shape  of  the  cross  section  was  assumed 
rectangular  for  all  stiffeners. 

3.2.  LATERALLY  LOADED  CYLINDRICAL  SHELL  STIFFENED  WITH  NON-UNIFORM  RINGS 

3.2.1.  APPROXIMATE  ANALYSIS 

Lateral  pressure  loading  is  of  less  practical  importance  than  hydrostatic  pressure  but  is  easier 
to  handle  analyticly.  Hence  it  is  considered  first.  ’Now  general  instability  under  lateral  pressure  oc¬ 
cur*,  for  unstiffened  cylindrical  shells  as  well  as  stiffened  ones,  with  one  wave  in  the  axial  direction, 
sec  [1]  and  [16].  Hydrostatic  pressure  is  a  case  of  combined  loading  (axial  and  lateral  pressures)  and 
yields  buckling  shapes  that  are  combinations  of  the  modes  occuring  under  the  axial  and  lateral  compo¬ 
nents.  This  interaction  was  analysed  for  stiffened  shells  in  [1]  for  a  large  range  of  shell  and  ring  ge¬ 
ometries.  It  was  found  there  that  for  short  and  thick  shells  (small  Z)  general  instability  occured  with 
many  waves  in  the  axial  direction.  From  these  arguments  one  can  qualitatively  conclude  that,  while  for 
lateral  pressure  rings  of  non-uniform  cross  section  varying  linearly  or  according  to  half  a  sine  wave, 
see  Fig.  2,  might  be  very  promising,  for  hydrostatic  pressure  the  same  distribution  of  rings  will  inevi¬ 
tably  lead  to  large  losses  of  efficiency.  This  point  i-.  later  discussed  in  detail.  For  cylinders  stiffen¬ 
ed  with  rings  only,  the  first  set  of  N  algebraic  equations  in  Eqs.  (11)  uncouple,  and  hence  the  order 
of  the  stability  determinant,  Eq.  (17),  is  reduced.  For  more  details,  as  well  as  explicit  expressions  of 
integrals  (16),  see  Appendix  D. 

To  get  an  idea  of  the  structural  efficiency  of  cos-uniform  stiffeners,  the  csss  of  rings  with  non- 
uniform  width  ,i  first  analysed.  A  closed  form  solution  is  possible  if  one  uses  the  first  terms  of  Eqs.  (9) 
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as  a  rough  approximation  to  the  buckling  shape  for  general  instability. 

For  a  width  of  ring  varying  according  to  f(x)  along  the  generator,  the  weight  saving,  relative 
to  a  uniformly  stiffened  shell  buckling  at  the 'Same  load,  is  given  by 

=  [l  +  (A2/.h)l/(l  +  (kA2/ah>] 

where 

(«•// 3)  (ir/i 3) 

k  *  (2  f  f(x)  sin2  /3xdx)/(  f  f(x)  dx)  (33) 

o  o 

and  A2  is  the  area  of  rings  for  the  uniformly  stiffened  cylinder. 

The  function  f-x)  which  minimizes  the  ratio  (largest  weight  saving)  was  found  to  be  Dirac's  delta 
function. 

f(x)  -  S[OrA0)-x]  (34) 

for  which  the  value  of  the  ratio  is 

h^./h^  =  [I  +  (A2/«h)]/[l  +  2(A2/ah)l  (35) 

With  A£/  ah  =  1 ,  which  represents  a  typical  geometry  in  aerospace  vehicles,  a  maximum  weight  sav¬ 
ings  of  33.3%  is  then  obtained. 

Variation  of  ring  width  according  to  Eq.  (34)  means  concentration  of  the  whole  stiffener  area  in 
one  bulkhead  at  ‘he  mid-length  of  the  shell,  see  Fig.  3.  This  is  the  direct  result  of  the  approximation 
es.Mmed  at  the  beginning  of  this  section.  The  deflecting  shape  for  general  instability  assumed  a 
half  sine  wave  (n  =  1),  hence  the  rings  tend  to  concentrate  in  the  weakest  part  of  the  shell  -  in  the 
middle.  The  n-2  buckling  mode,  however,  was  not  taken  in  account,  the  ugh  in  the  case  presented  in 
Fig.  3  it  is  the  one  most  prone  to  occur.  The  shell  will,  therefore,  buckle  locally  at  a  pressure  much 
below  that  causing  general  instability  of  the  uniformly  stiffened  cylinder.  The  efficiency  o i  con-uniform 
stiffeners  is  therefore  made  up  of  two  opposing  contributions.  A  primary  positive  contribution  resulting 
from  the  strengthening  of  the  weak  mid-length  of  the  shell,  and  a  secondary  negative  contribution  due 
to  the  weakening  of  the  shell  close  to  tbs  edges  that  may  thee  buckle  prematurely.  This  premature 
buckling  mode  will  henceforth  be  called  a  "many  wave"  deflection  shape  in  the  axial  direction. 
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NON-UNIFORM  HEIGHT 

obtain  a  better  feeling  for  the  influence  cf  the  opposing  contributions  on  structural 
efficiency,  different  types  of  height  variation  were  analysed.  In  Fig.  4  the  non-dimensional  ratio 
c2/c2  that  describes  the  variation  of  the  height  of  the  rings  along  a  generator,  is  drawn  versus  x 
for  different  values  of  the  parameter  y.  y  the  "variation  ratio"  is  the  ratio  between  the  weight  of 
the  uniform  part  of  the  ring  to  their  total  weight  and  its  values  are  restricted  by  the  inequality 

0  *  y  <  1  (36) 

y  =  1  represents  the  case  of  uniform  stiffeners,  while  y  «  0  the  one  where  the  entire  height  of  the 
ring  varies  *  m wsa^aUy  with  x.  The  smaller  y  the  more  of  the  stiffeners  weight  (which  is  kept  con¬ 
stant  for  all  the  values  of  y)  is  near  the  mid-length  of  the  shell. 

In  Fig.  5  the  radial  displacement  at  buckling  is  given  for  a  short  shell  and  different  values  of 
y.  The  parameters  e2/h  and  A^/ah  are  those  of  the  equivalent  uniform  rings  (weight  equivalence). 
This  equivalency  will  he  used  later  on  for  computing  the  efficiency  of  non-uniform  stiffeners.  For 
y  =  1  the  deflection  is  a  half  sine  wave;  a  decrease  in  y  distorts  the  buckling  shape  so  as  to  flatten 
its  extremum  at  x  -  (L/2ft).  It  is  worth  mentioning  that,  except  slight  corrections  in  shape,  the  curve 
still  resembles  the  basic  half  sine  wave  in  the  axial  direction.  This  supports  the  intuitive  argument 
given  in  Section  2-3.  As  y  is  further  decreased,  however,  a  change  ip.  behaviour  is  observed.  For 
y  ■  0.25  two  completely  different  deflection  curves  are  drawn.  One  represents  the  normal  shape  with 
one  buckle  along  the  generator,  whereas  the  other  represents  a  "multi-wave"  behavior  with  relatively 
small  displacements  at  the  middle  of  the  shell.  During  the  process  of  minimization  of  the  lowest  eigen¬ 
value  of  (17)  relative  to  t,  a  double  minimum  curve  is  found  (see  for  example  Fig.  9) .  Each  of  the 
extremus  represents  a  different  buckling  mode,  and  the  one  yielding  the  lowest  buckling  pressure  is 
physically  significant  In  Fig.  5  both  deflection  curves  were  drswn  for  y  «  0.25  because  they  both 
occured  at  almost  the  same  value  of  the  externei  pressure.  (The  difference  in  A  was  of  the  order 

P 

of  0.9%).  For  y«0  the  n/1  buckling  mode  is  dominant. 

In  Figs.  <3  and  7  the  radial  displacements  are  given  for  two  other  shells  of  different  geometry. 
The  longer  rhe  shell  the  more  pronounced  is  the  mode  change.  The  term  "edge  buckling"  can  be  ap- 
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In  order  to 


18 


plied  to  the  "multi-wave"  buckling  pattern  in  the  aenae  that  very  rapid  changes  in  the  displacements 
take  place  in  the  neighbourhood  of  the  edges.  The  magnitude  of  the  displacements  is  zero  for  more  than 
70%  of  the  length  of  the  shell.  ‘In  Fig.  7  two  curves  are  again  drawn  for  y  ■=  0.25  for  the  same  reason 
as  before  in  Fig.  5. 

In  Table  1  the  computed  eigen-vectora  are  given  for  the  deflection  shapes  of  Fig.  7.  For  y  «  0 
and  0.25  the  two  columns,  represent  two  possible  buckling  shapes.  The  first  is  the  so-called  "n  =  1" 
mode  of  instability.  This  name  is  fully  justified  since  Cj  is  much  larger  than  the  other  coefficients. 

The  second  column  represents  the  n  1  buckling  shape.  It  is  characterised  by : 

a.  poor  convergence  of  the  stability  determinant,  Eq*  (17); 

b.  a  large  number  of  circumferential  waves  (see  table  1) ; 

c.  none  of  the  terms  of  the  displacement  series,  Eqs.  9,  dominates. 

For  the  particular  case  y»  0  the  terms  C9  ,  Cjj  and  CJ3  are  of  similar  magnitude. 

In  Fig.  8  the  deflection  curves  for  the  same  geometries  of  shell  and  rings  as  in  Fig.  6  are  given 
for  outside  rings  (e2/h  <  0).  There  is  a  noticeable  influence  cf  the  eccentricity  on  the  efficiency  of 
the  non-uniform  variation  of  stiffeners.  A  comparison  of  Figs*  6  and  8  shows  that  for  outside  rings  the 
"multi-wave"  buckling  shape  is  postponed  ts  lower  values  of  y.  As  will  be  seen  later,  this  increases 
the  gain  in  load. 

Figs.  5  to  8  clearly  show  that  from  a  certain  value  of  y  (the  actual  value  depends  on  the  shell 
geometry)  down,  an  "edge  buckling"  pattern  dominates,  leaving  the  middle  part  of  the  shell  relatively 
underformed.  Numerical  results  for  this  "edge  buckling"  mode  show  large  losses  in  efficiency.  The 
buckling  load  amounts  to  less  than  50%  of  that  obtained  with  uniform  stiffeners.  The  physical  argumeots 
behind  those  losses  are  similar  ia  nature  to  those  presented  in  the  analysis  of  Eq,  (34).  in  other  words, 
the  oa*,rstiffening  of  the  middle  part  of  the  shell  leaves  the  neighbourhood  of  the  edges  relatively  weak 
and  therefore  very  prone  to  buckling. 

In  Fig.  10  tbe  gain  in  load  pM/pu  is  drawn  versus  y  for  different  shell  geometries.  Most  of 
the  curves  have  a  discontinuity  in  slope  at  certain  values  of  y  that  may  be  denoted  yfc.  'For  y  >  yfc 
the  deflections  during  buckling  have  one  basic  wave  along  the  generator.  Decrease  in  y  from  y  =  1 
*°  ftr.  rc*ull,‘  *n  8  monotonous  increase  of  the  gain  in  ioad.  fo  the  left  of  the  discontinuity  the  "edge 
buckling"  pattern  dominates,  leading  to  losses  in  efficiency. 


i 
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For  y  -  on  j  obtains  ihe  highest  gains  in  load  for  a  certain  shell  geometry*  Its  value  varies  with 

Z.  For  very  short  shells,  for  example*  no  *multi*wave*  pattern  takes  place  and  therefore,  *  0  is 
the  optimum  variation  ratio* 

In  Fig.  11  the  optimum  gain  in  load,  obtained  at  y  -  y  (see  Fig.  10);  is  plotted  versus  the 
shell  geometry  parameter  Z.  For  very  low  values  of  Z,  Z  <  30.  is  zero.  In  the  intermediate 
range  30  <  Z  <  100  very  rspid  changes  occur  in  the  value  of  y#  csuaing  the  changes  in  curvature 
that  appear  in  Fig.  11.  Further  away  from  the  perturbed  region  the  rise  of  the  curve  is  monotonous.  It 
should  be  noted  that  for  long  shells,  the  mode  transition  is  very  sharp,  sea  Fig.  10,  and  causes  very 
large  losses  in  efficiency.  A  small  msrgia  of  safaty  is.  therefore,  advisable  and  variatioa  ratios  slight* 
ly  larger  than  y*  should  be  employed.  The  larger  gains  in  load  obtained  with  outside  rings  are  some¬ 
what  overshadowed  by  the  superiority  of  uniform  inside  rings  over  outside  oaee.  see  [11 . 

In  most  problems  of  practical  interest  the  geometry  of  the  ahull  as  well  as  the  applied  loads  are 
given.  I  he  designers'  job  is  to  find  ths  configuration  of  stiffeners  yielding  minimum  weight.  Care 
should  be  taken  about  physical  end  practical  constraints  such  as  yielding  ot  material,  maximum  permis¬ 
sible  ring  or  stringer  height,  etc. 

The  procedures  deacribed  in  Chapter  2  far  sefvieg  differ  serial  aq**rieas  (5)  at*  efcea  known  as 
the  "direct"  problem:  Computation  of  the  load  causing  general  instability  for  given  stiffener  and  shell 
geometry.  The  "indirect*  problem  of  finding  the  ring  geometry  for  e  given  load  is  in  general  more  dif¬ 
ficult.  One  assumes  different  values  for  the  ring  geometry  end  computes  the  corresponding  values  of 
the  buckling  load  with  one  of  the  methods  of  Chapter  2.  When  these  values  ere  close  enough  to  the  load 
given  in  the  data,  a  regular  falsi  method  is  used  to  increase  the  convergence  of  solution. 

In  Fig.  11  the  weight  saving  equivalent  of  the  gain  load,  given  by  the  full  curves,  is  plotted 
versus  Z.  its  magnitude  is  never  larger  then  10-15%  .  It  should  be  carefully  noted  that  the  present 
analysis  is  not  an  optimization  procedure  in  the  sense  of  [17]  end  [18]  because : 

a.  only  general  instability  of  the  composite  structure  is  analysed, 

b.  the  weight  minimization  ie  on  stiffener  configuration  only,  whereas  the  thickness  of  shell  is  kept 
constant 

in  rig.  11  the  weight  of  the  rings  t*  only  1/3  of  the  tolsl  weight  (A2/a  =  h/2).  Therefore,  saving 
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15%  in  the  total  weight  of  the  shell,  by  using  non-uniform  stiffeners,  is  equivalent  to  a  40%  saving  in 
the  weight  of  the  rings. 

The  structural  efficiency  depends  to  a  larg  extent  on  the  area  ratio  of  the  rings  (A2/ ah),  and 
the  larger  the  area  ratio  the  larger  are  the  weight  savings.  Typical  computations  have  shown  that  for 
(A2/ah)  =  1  they  amount  to  20—25%  of  the  original  weight  of  the  uniformly  stiffened  shell. 

In  Fig.  12  and  Table  2  the  gain  in  load  is  plotted  versus  Z  for  different  values  of  y.  The 
first  two  curves  for  y  =  0.75  and  9.5  show  a  continuous,  monotonous  rise  of  efficiency  with  Z.  The 
radial  displacements  at  buckling  have  one  lobe  in  the  axial  direction.  For  y  =  0.25  the  curve  for  in¬ 
side  rings  has  two  discontinuities  in  slope.  'In  the  intermediate  range  of  Z  ,  the  "edge  buckling"  pat¬ 
tern  dominates  and  causes  leases  in  efficiency.  It  should  be  mentioned  that  the  curve  for  outside  rings 
and  the  same  value  of  y  is  continuous.  The  curves  for  y  *  0  are  mainly  of  academic  interest,  since 
they  emphasize  the  consequences  of  understiffened  edges  on  general  instability. 

In  Fig.  13  and  14  the  influence  of  ring  cross-sectional  area  A2  and  magnitude  of  eccentricity 
(e2/h)  on  the  gain  in  load  is  investigated.  Comparing  respective  curves  in  Fig.  12  and  13  shows 
small  differences  in  gain.  On  the  other  hand,  the  curves  in  Fig.  14  are  of  somewhat  different  nature, 
fhe  "edge  buckling"  pattern  appears  already  for  y  =  0.5,  while  for  y  =  0.25  its  zone  of  influence 
extends  to  a  larger  range  of  Z.  This  is  not  surprising  since  increasing  the  equivalent  height  of  rings 
from  (e2/h)  =  5  to  (r„/h)=10  results  in  overstiffening  the  middle  part  of  the  shell.  A  clearer  picture 
is  obtained  in  Fig.  15  where  the  gain  is  plotted  for  different  values  of  (e2/h)  and  a  constant  y.  For 
(e2/h)  =  3  the  curve  is  continuous  with  relatively  high  gains  for  all  the  spectrum  of  Z.  For  the  other 
values  of  (e2/h)  the  efficiency  gradually  decreases  and  may  even  change  i’s  sign.  While  Fig.  15  anal¬ 
yses  the  influence  of  (e2A)  on  the  ’edge  buckling"  pattern.  Fig.  16  does  the  same  for  a  pattern  with 
one  buckle  along  the  x  axis.  For  long  shells  (large  Z)  Fig.  16  shows  larger  gains  in  load  the  larger 
the  eccentricity.  In  the  left  part  of  the  figure  (small  Z)  the  behavior  of  the  curves  is  simitar  that  in 

Fig.  15. 

3.2.3.  RINGS  OF  NON-UNIFORM  WIDTH 

When  comparing  a  height  and  width  variation  of  rings,  the  superiority  of  the  first  is  almost 
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obvious,  since  the  stiffening  material  is  then  used  in  a  more  efficient  way*  'In  many  cases,  however, 
practical  constraints  rule  out  the  use  of  height  variation,  leaving  only  width  variation  as  a  possible 
method  of  optimization.  The  explicit  formulation  of  integrals  (16)  for  a  sinusoidal  width  variation  is 
given  in  the  third  section  of  Appendix  D.  The  absence  of  stringers  simplifies  Lqs.  (12)  yielding 
2N  x  2N  stability  determinant.  The  algebraic  details  are  similar  to  those  given  for  height  variation 
(see  Appendix  C). 

Numerical  results  show  that  in  the  present  case  a  variation  ratio  y  =  0  yields,  at  buckling,  a 
deflection  shape  with  one  lobe  in  the  axial  direction.  Computations  were  carried  out  for  a  large  spec¬ 
trum  of  shell  and  ring  geometries? 

For  larger  values  of  y  the  same  mode  persists  leading  to  a  gradual  decrease  in  efficiency. 
Hence  only  the  optimum  y  is  used  in  the  computations.  It  may  be  recalled  that  in  the  case  of  height 
variation  of  rings,  y  =  0  yielded  very  low  efficienci  es.  If  one  aims  at  a  physical  explanation  of  the 
difference  between  both  variations,  one  finds  that  a  non-uniform  width  understiffens  the  edges  in  a 
milder  way  postponing  the  appearance  of  the  "edge  buckling*  pattern.  In  the  mathematical  formulation, 
the  difference  between  width  or  height  variation  according  to  sin  /3x  lies  in  the  fact  that  the  bending 
stiffness  of  the  ring  changes  according  to  sin  /3x  and  sin3  /3x  respectively.  The  intensity  of  the 
non-uniform  distribution  at  the  edge  of  the  shell,  given  by  the  slope  of  the  sin  )3x  and  sin3  /9x  curves 
is,  therefore,  finite  for  width  variation  and  zero  for  height  variation.  Since  this  slope  is  of  critical  im¬ 
portance  for  the  "edge  buckling"  pattern,  the  difference  in  behavior  is  apparent.  In  Fig.  17  and  Table 
2  the  gain  in  load  and  weight  saving  are  plotted  versus  Z  for  inside  and  outside  rings.  Comparison 
of  the  results  with  the  optimum  gains  obtained  for  height  variation,  see  Fig.  11,  confirms  the  intuitive 
conclusions  reached  at  the  beginning  of  the  section.  In  Fig,  17  the  maximum  gains  in  load  are  of  the 
order  of  15  —  20%,  much  lower  than  those  for  the  height  variation.  The  weight  savings  on  the  other  hand 
are  of  the  same  order  of  magnitude.  This  is  not  surprising  since  the  weight  optimization  was  performed 
in  each  case  with  different  constraints.  For  rings  of  non-uniform  width,  the  efficiency  of  stiffeners  was 
compared  with  that  of  uniform  rings  of  the  same  height  and  equivalent  width.  For  height  variation  the 
equivalency  was  related  to  the  height,  while  the  width  was  kept  constant. 
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3.3.  AXIALLY  LOADED  CYLINDRICAL  SHELL  STIFFENED  #Y  NON-UNIFORM  STRINGERS 

Recent  optimization  studies  on  axially  loaded  cylindrical  shells  [10]  and  [18]  have  shown  that 
a  combination  of  rings  and  stringers  yield  the  highest  structural  efficiency*  In  the  present  analysis 
only  the  influence  of  non-uniform  stringers  is  studied  in  detail-  That  of  rings  alone, is  of  no  interest, 
since  alone  they  are  very  inefficient  for  this  type  of  loading*  An  optimization  with  both  non-uniforr 
stringers  and  rings  could  yield  further  gains,  but  is  not  considered  in  this  report.  Sinusoidal  and  linear 
height  variations  (sea  Fig.  2)  were  analysed  in  order  to  study  the  influence  of  the  functional  variation 
on  structural  efficiency.  ‘A  width  variation  of  stringers  is  of  little  interest  since  it  yields  gains  in  load 
of  only  few  percents.  In  the  case  of  stringer-stiffened  shells  every  term  of  the  displacement  series, 

Eqs.  5',  is  a  possible  solution  of  the  secorn.  equation,  Eqs.  (5).  By  relatively  simple  manip¬ 

ulations,  the  order  of  the  matrix,  Eq.  (17),  can  be  reduced  to  2N .  Further  details,  as  well  as  the  ex¬ 
plicit  expressions  of  integrals  (16)  for  the  two  types  of  variation,  are  given  in  Appendix  E. 

In  Fig.  18  and  Table  3  the  influence  of  shell  geometry  on  the  structural  efficiency  of  non-uniform 
stringers  is  studied  for  different  values  of  y.  For  y  -  0.75  and  0.3  the  curves  drop  monotonously 
to  asymptotic  value  P^/p,^  =  1*  For  long  shells  (Z  >  5000)  stiffened  by  inside  stringers,  the  gains 
in  load  do  r.ot  justify  the  increase  in  the  cost  of  production  due  to  non-uniform  stringers.  It  should  be 
mentioned,  that  in  the  same  range  of  Z ,  even  uniform  stiffeners  are  less  efficient  then  equivalent 
thickening  of  shell,  nee  [3]  •  Comparing  the  gains  in  i'.oad  for  ring-stiffened  shells,  Fig.  12,  and  stringer- 
stiffened  ones.  Fig.  18,  one  observes,  in  addition  to  the  difference  in  magnitude,  that  the  behavior  with 
increase  in  Z  is  reversed. 

For  y  «  0.25,  and  even  more  so  for  y»  0  (see  Table  2),  an  "edge  buckling"  pattern,  simile  - 
to  the  one  occuring  in  ting  stiffened  shells,  dominates,  causing  a  very  rapid  decrease  in  structural  ef¬ 
ficiency.  The  reduction  is  even  more  pronounced  than  that  occuring  in  Fig.  12,  because  of  the  different 
buckling  behaviors  of  isotropic  shells  under  lateral  and  axial  loadings.  For  large  values  of  Z  (long 
shells),  the  single  "half  wave"  buckling  pattern  is  once  again  dominant  (even  for  y«0),  yielding  higher 
gains  for  lower  values  of  y. 

For  small  Values  of  Z  (Z  <  100),  cylindrical  shells  stiffened  with  outside  stringers  buckle  in 
an  axisymmetric  mode,  see  also  [3],  causing  a  large  scatter  of  ti  :  computed  poiuts.  The  curves  in  Fig. 
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58  smooth  out  he  effect  of  this  scatter,  In  Fig,  18  the  influence  of  eccentricity  on  structural  efficiency 
is  also  investigated.  For  values  of  Z  smaller  than  1000,  the  gains  in  load  are  higher  for  inside 
stringers.  An  inversion  occurs  at  larger  Z  values,  yielding  noticeable  differences  between  the  effi* 
ciency  of  both  configurations*  ‘For  a  typically  long  shell,  (L/R)  =2,(R/h)  =  1000,  the  gain  in  load  is 
increased  from  2%  to  15%  by  changing  the  position  of  the  stiffeners*  It  should  also  be  mentioned  that 
recent  experimental  and  theoretical  investigations,  [3]  to  [61,  emphasize  the  superiority  of  uniform  out¬ 
side  stringers  over  inside  ones*  ‘In  Fig.  19  the  influence  of  the  magnitudd  of  eccentricity  on  the  gain 
in  Lad  is  investigated*  Since  a  variation  ratio  y«0*3  is  used,  "edge  buckling*  is  eliminated  as  a 
possible  mode  of  instability*  Increase  in  the  magnitude  of  eccentricity  for  a  constant  string  sr  area, 
yields  more  slender  stiffeners,  whose  height  increases  in  inverse  proportion  to  the  decrease  of  the  width* 
As  a  result,  the  moment  of  inertia  of  stiffener  cross  section  varies  according  to : 

<lu/bh3)«(ei/h)2 

The  gains  in  load,  which  are  in  direct  proportion  to  the  bending  stiffeness  of  the  stringers,  grow 
larger,  the  larger  the  eccentricity*  However,  practical  considerations  prescribe  upper  limits  to  the  mag¬ 
nitude  of  e./h*  Some  of  them  are: 

a)  local  buckling  of  stiffener  as  a  simply  supported  -  free  plate, 

b)  the  thickaess  of  commercially  available  sheets, 

c)  reduction  of  the  width  of  the  stringer  dj  (see  Fig.  1)  beyond  a  certain  limit,  reduces  the 
"effective  width"  of  the  sheet  between  the  stiffeners*  The  present  analysis  does  not  account  for  the 
resulting  decrease  in  the  total  stiffeness  of  the  shell. 

In  Fig*  20  the  influence  ef  y  on  the  gain  in  load  is  investigated  for  two  types  of  non-uniform 
height  variation.  In  the  neighbourhood  y  *■  1 ,  the  two  variations  yield  similar  results*  With  decrease 
in  y  the  gain  in  load  increases  and  goes  through  a  maximum,  which  represents  a  general  characteristic 
behavior  of  non-uniform  stiffeners.  In  the  case  of  rings,  this  maximum  appeared  as  s  discontinuity  in 
slope,  see  Fig.  10,  This  discontinuity  is  caused  by  transition  from  an  "edge  buckling"  pai'em  to  a 
single  half  wave- pattern  in  the  axial  direction*  'In  the  present  case  of  stringers,  however,  the  transition 
is  continuous*  For  sinusoidal  height  variation  the  optimum  value  of  y  is  lower  than  for  the  lineer  one, 
yielding  therefore  higher  gains  in  load.  Physically,  this  means  that  the  linear  variation  is  more  prone 
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to  "edge  buckling*  and  "needs*  larger  parts  of  the  stringer  to  be  uniform  along  the  shell.  The  optimum 
y  at  which  the  maximum  gains  are  obtained  varies  with  the  shell  geometry.  For  long  shells  (large  Z) 
it  has  o  tendency  to  move  to  the  left  of  Fig.  20  until  a  value  of  Z  is  reached  (Z.«  10  )  for  which 

3 

y  -  0  is  the  best  configuration.  In  the  "short  shell*  range  Z  <  10  higher  values  of  y  are  needed 
to  increase  the  stiffness  in  the  axial  direction. 

In  Fig.  21  the  two  types  of  stiffener  variation  are  compared  for  different  shell  geometries.  While 
for  y  =  0.75  the  linear  variation  is  superior,  the  opposite  holds  for  y  =  0.5  and  0.25. 

For  axially  compressed  isotropic  thin  cylindrical  shells  the  experimental  buckling  stresses  are 
much  below  those  computed  by  small  deflection  theory.  Recent  experimental  results  [4]  [6]  and  [19] 
have,  however,  shown  that  heavily  stiffened  shells  can  be  adequately  analysed  by  linear  theory.  With 
closely  spaced  and  relatively  large  stiffeners,  experimental  reaults  fall  within  a  few  percent  of  the 
theoretically  computed  loads.  For  small  values  of  y  the  edges  of  the  shell  are  relatively  undei-stiffen- 
ed,  and  hence  the  applicability  of  linear  theory  is  in  doubt.  'Therefore,  the  variation  ratio  y  =  0.5 
seems  preferable,  even  though  for  long  shells  y  *>  0.25  yields  higher  structural  efficiencies  than 
y  =  0.5,  (see  Fig.  18). 

3.4  COMPARISON  BETWEEN  THE  GALERKIN  METHOD  AND  THE  METHOD  OF  "CORRECTING 

COEFFICIENTS" 

The  nuraericui  results  of  the  preceding  sections  were  obtained  by  application  of  a  Galerkin  pro¬ 
cedure  to  Eqs.  (5).  The  stringer-stiffened  axially  compressed  cylindrical  shell  is  solved  again  by  the 
method  of  "Correcting  Coefficients",  the  details  of  which  are  given  in  the  theoretical  analysis.  'Tie 
results  are  compared  to  the  more  accurate  ones  obtained  in  the  Inst  section. 

The  expressions  of  Eqs.  (31)  were  used  for  the  correcting  coefficients,  with  the  implied  assump¬ 
tion  that  the  radial  displacements  at  buckling  have  one  half  wave  in  the  axial  direction.  This  assump¬ 
tion  is  incorrect  for  low  values  of  y,  where  the  "edge  buckling*  pattern  dominates. 

In  Fig.  22  and  Tabic  4,  the  ratio  P^/P^ ,  as  obtained  by  the  two  different  methods,  is  plot¬ 
ted  versus  Z  for  different  values  of  y.  For  y  »  0.75  the  maximum  difference  amounts  to  less  then 
I  ;  but  for  y  =  0.5  it  is  6%. 

In  general,  the  smaller  y  the  larger  the  differences.  Typical  cases  computed  with  y  »  0 
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showed  differences  of  more  than  100%.  The  diecrepancies  result  from  a  poor  choice  of  the  correcting 
coefficients.  A  better  approximation  (which  was  not  used  in  the  present  work)  employs  Eqo.  (32)  for 
computation  of  the  coefficients.  Mathematically  speaking,  the  method  of  "correcting  coefficients"  has 
no  rigorous  proof  of  convergence.  It  is  based  on  physical  and  intuitive  arguments.  In  [15]  where  the 
method  was  used  for  the  first  time,  the  results  were  not  compared  to  those  of  a  more  rigorous  solution. 

An  approximate  evaluation  of  their  accuracy  is  obtained  from  the  results  of  the  present  solution. 

The  main  features  of  the  method  are : 

a)  reduction  of  the  order  of  stability  determinant  of  Eq.  (17). 

b)  quicker  convergence  for  the  same  number  of  terms  in  series  (27)  than  the  more  accurate  solution. 

c)  the  values  of  the  loads  it  yields  are  non*conservative, 

d)  its  accuracy  depends  to  a  large  extent  on  the  choice  of  the  correcting  coefficients. 

3.3,  CYLINDRICAL  SHELL  UNDER  HYDROSTATIC  PRESSURE  AND  STIFFENED  WITH  UNIFORM 
STRINGERS  AND  RINGS  OF  NON-UNIFORM  HEIGHT. 

Uniformly  ring-stiffened  shells  under  hydrostatic  pressure  may  buckle  in  one  of  two  possible 
shapes,  see  [1].  For  short  shells.  Z  <  1000,  and  relatively  heavy  stiffeners,  the  n  ^  1  bucxiing 
shape  is  dominant  reducing  the  efficiency  of  rings  drastically.  In  [10]  stringers  were  combined  with 
rings  to  improve  the  structural  efficiency.  It  was  rhown  there  that  longitudinal  stiffeners  of  cross  sec¬ 
tional  area  amounting  to  only  5  —  10%  of  that  of  the  rings,  increase  the  buckling  load  of  the  shell  by 
50-70%.  Further  increase  in  the  stringer  area  (with  corresponding  smaller  rings)  reduced  the  buckling 
pressure,  because  once  the  o  =  1  mode  of  instability  is  dominant,  rings  are  the  best  stiffeners. 

The  present  analysis  combines  uniform  stringers  with  rings  of  non-uniform  height.  The  uniform¬ 
ity  of  the  stringers  simplifies  the  expressions  of  the  integrals  (16),  reducing  thereby  the  order  of  stabil- 
it/  matrix,  Eq.  (17).  Further  detail#  are  given  in  Appendix  F . 

Fig.  23  and  Table  5  analyse  the  influence  of  sinusoidal  height  variation  of  rings  on  the  critical 
hydrostatic  pressure.  For  values  of  Z  below  1000,  negative  gsins  are  obtained.  These  losses  in  ef¬ 
ficiency  may  be  attributed  to  the  two  main  causes: 

s)  For  snort  shells,  Z  <  lOOo  Jformly  stiffened  shells  buckle  in  an  n  p  1  deflection  shape. 
With  non-uniform  rings,  the  overstiifening  of  the  middle  of  the  shell,  at  the  expeoee  of  the  edges,  causes, 
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therefore,  large  losses  of  efficiency  in  this  range  of  Z . 

b)  For  lateral  pressure,  an  "edge  buckling*  pattern  dominates  for  certain  values  of  y  and  shell 
geometry.  A  similar  effect  occurs  in  the  case  of  hydrostatic  pressure  loading,  and-reduces  the  efficien¬ 
cy  f  f  non  uniform  stiffening. 

These  effects  will  henceforth  be  referred  to  as  the  *first"  and  "second"  effect. 

For  short  shells  with  Z  below  100,  Fig.  23,  the  first  effect  dominates,  and  there  is  practically 
no  g .in  in  load.  For  longer  shells  the  curves  go  through  a  minimum  whose  magnitude  depends  primarily 
on  the  second  effect*  The  lower  the  y  the  larger  the  inefficiency  of  non-uniformed  stiffening.  '  Note 
that  the  curves  to  the  right  of  the  discontinuities  in  slope  in  Fig.  23,  represent  n  =  1  buckling  shapes. 
The  magnitude  of  the  gain  is  then  similar  to  that  obtained  for  lateral  loading,  Fig.  12 .  In  order  to  ex¬ 
plore  possible  improvements  in  the  structural  efficiency  of  shells  (for  stability),  different  combinations 
of  stringers  and  rings  (of  constant  weight)  were  analysed.  The  results  are  presented  in  Fig.  24  for  a 
variation  ratio  y  *0.5.  '  Results  for  other  values  of  y  are  given  in  Table  5*  Each  curve  in  Fig.  24 
has  two  discontinuities  in  slope.  ‘The  intermediate  range  of  Z,  for  which  the  n  £  1  buckling  shape 
dominates,  decreases  gradually  with  increasing  stringer  area.  ‘For  large  values  of  Z  all  curves  merge 
into  one. 

If  one  aims  at  structural  optimization,  one  must  remember  that  increase  of  stringer  area  beyond 
a  certain  limit  is  detrimental  to  the  efficiency  of  the  stiffeners.  The  values  of  (A^)^  for  different 
stiffener  configurations  are  compared  in  Table  5  for  a  typical  long  shell  (L/R=2,  R/b=1000).  A  re¬ 
duction  of  25%  in  buckling  load  is  found.  For  shorter  shells  an  increase  appears  instead,  but  it  gradual¬ 
ly  changes  its  sign  to  a  reduction  with  increase  in  stringer  area.  Hence  an  optimum  value  of  stringer 
area  exists  for  every  shell  and  ring  geometry.  Comparison  of  Fig.  25  with  Fig.  24  brings  out  the  in¬ 
fluence  of  eccentricity  of  rings.  The  comparison  shows  that  with  outside  rings,  stringers  of  smaller 
cross-sectional  area  are  needed.  A  similar  behavior  was  already  noticed  in  Fig.  12,  where  outside  rings 
were  also  found  to  be  better. 

For  hydrostatic  pressure  loading,  both  the  symmetric  and  anti-symmetric  components  of  the  sta¬ 
bility  matrix,  Eq.  17,  have  to  be  checked  carefully.  Axisymmetric  buckling  may  occur  for  short  shells 
with  outside  rings  (see  Table  5) .  Finally,  one  should  remember  to  check  the  curve  of  load  versus  num¬ 
ber  of  circumferential  waves  (t)  in  each  case  for  a  double  minimum  (see  Fig.  9). 
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APPENDIX  -  A 

THE  GEOMETRICAL  PARAMETERS  OF  THE  STIFFENERS 


Explicit  expressions  foe  the  geometrical  parameters  of  the  stiffeners  are  given  below,  first  for 
a  general  stiffener  cross-section  and  then  for  stiffeners  of  rectangular  cross-section.  The  distance 
between  stiffeners  is  assumed  to  be  uniform. 


P  arametcr 

General  Definition 

Rectangular  Cross  Section 

MjW 

0  -  t>2)  Aj  (x)/bh 

(1  -  t/2)[c,(x)/h]  (dj(x)/b] 

(1  -  v”)  A2  (x)/ah 

(1  -  v1)  [  c2  (x)Al  W2  (x)/a] 

/t*(x)  et(x)/R 

ftj(x)  l[cj(x)  +  h]/2B  I  < 

y*(x) 

(x)  e2  (x)/R 

#*J(x)  l[c2(x)  +h]/2R}r 

<TT<x) 

12(R/h)2vJ(x) 

C2M 

12(R/h)Jy*U) 

12(1  -u2)  IQ1  (x)/bh* 

12  b  4  b 

,?0  2(x) 

12(1  -  »/2)  I02(x)/«b3 

<1  (x> 

12(1  -  v2)  GI^j/Eb!,3 

2 

[2/(1  +  v)]  M*(x)  [(dj  (x)/b)  (b/h)] 

<2(x) 

12(1  -  i/s)  GI^/Ebb* 

{2/(l  +  ^)]>u*(x)[(d2(xVe)(«A>]i 

where 

for  inside  stiffeners 
for  outside  stiffeners 

Note  that  in  the  last  two  expressions,  the  torsion  constant  of  a  rectangular  cross-section  is  computed 
with  the  assumption  of  c/d  -  5  ,  Hence  I,  ,  (x>  =  {  c ,  (x)  dj  (x)]/2. 
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APPENDIX  -  0 

THE  INTEGRALS  I((n,m)  AND  K^n.m)  FOR  NON-UNIFORM  STIFFENERS  WITH  A  DISCONTINUITY 
IN  SLOPE  (im  Fig.  24). 

Evaluation  of  integral*  (16)  by  bisection  of  tbe  limits  of  integration  is  incorrect.  The  separa¬ 
tion  must  be  carried  out  from  Eq.  (5),  by  an  analysis  of  two  sets  of  differential  equations  for  the  differ¬ 
ent  domains. 

Another  approach  to  the  problem  uses  tbe  expression  of  the  virtual  wo/k  instead  of  Eq.  (1). 


2 n  (L/R) 


[  l  *  Vv  +  - w./RH&g  - (lyRKs^)  + 

>(M^/R)(S^)-(M^/R)(^)  +  Ns0«r;  +  N^05^]R^xd^  -  0 


(B-l) 


It  should  be  noted  that  neither  the  forces  nor  the  moments  in  Eq.  (B.  1)  are  differentiated  with 
respect  to  x . 

Now  t '  and  arc  non-linear  strains  defined  by 

<  -  («?,/»  <i  -  (8-2) 

Snbstitution  of  the  expression  for  the  strains  nod  curvatures,  Eqs.  (41  sad  (B.2),  into  (B.l) 


yields 


in  WR 


-(M^/RHSw)^  i-(Mx^Tl)(Sw)  ^  +  Ni0  w  M  (8w)  t  +  N^w  ^  (8w)  ^](R2dxd^)  -  0 


(B  -3) 


Integration  by  parts  of  Eq.  (B.3)  yields  the  equilibrium  equations,  Eqs.  (1),  and  the  boundary 
conditions  Eqs.  (?). 

Displacements  and  permissible  variations  are  assumed  as  in  Eqs.  (9)  snd  (10).  The  differentia¬ 
tion  and  integration  performed  is  shown  on  a  typical  term  of  Eq.  (B.3). 
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(B-y) 
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The  expression  for  I,  (n,m)  can  he  worked  out  in  a  similar  manner  as 
«+«  <V2R) 

I2(n,  m)  ■  (2/3/rr)  [1  +  (-1)  ]  /  %  *(x)  sin  n/3x  sin  m/9x  dx  (B-10) 


The  other  integrals  of  Eq.  (16)  are  evaluated  from  the  variational  expression,  Eq.  (B.  3),  in  a 


similar  manner. 


i 
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APPENDIX  -  C 

SOLUTION  OF  STABILITY  DETERMINANT. 

Application  of  Gauss's  algorithm  to  the  system  of  equations  (17)  yields 


where 

N 

* 

The  last  N 


is  the  symbol  for  aa  upper  triangular  matrix 
is  a  diagonal  matrix 
is  an  eigenvalue 
equations  uncouple,  yielding 


101 


(C-l) 


[M]1C  )  -  tffNJfC  1  (C-2) 

■  B 

premultiplication  by  M  yields 

(1/^)1CJ  -  (M-1]tNJ|C,l  -  [ffllfCJ  (C-3) 

[M]  is  a  non*sMigiitar  matrix,  as  otherwise  ^  »  0  is  an  eigenvalue  -  0  means  xero 
buckling  load). 

M  is  not  symmetric.  It  is,  however,  of  the  type  called  'symmetrisable",  since  the  original 
matrices  [M]  and  fNJ  are  symmetric  and  fNJ  is  positive  definite,  see  [20] .  The  eigenvalues  of 
M  are  therefore  all  real.  The  lowest  one  is  obtained  by  the  basic  technique  of  matrix  iteration.  The 
convergence  is  accelerated  by  exponential  extrapolation. 

It  should  be  mentioned  that  the  proposed  method  for  solving  the  stability  determinant  (17)  is 
quick  and  precise.  For  a  digital  computer  working  with  simple  precision,  6  digits  accuracy  was  ob* 
tained  for  matrices  of  order  ICC. 


32 


APPEND IX  -  D 

NON-UNIFORM  RING  STIFFENERS 

D.l  REDUCTION  OF  THE  ORDER  OF  STABILITY  DETERMINANT 
If  only  ring*  are  used  qs  stiffeners. 

H*M  -  q*,(x)  «•  X*(x)  «  r/t*(x)  =  0  (D.l) 

The  torsional  stiffness  of  the  rings  is  neglected  here,  and  therefore  (see  Eqs.  (16)) 

I,(n,aj)  =0  j  =  1.  2,  3.  4  (D.  2) 

The  matrices  [FA  ][FB1  and  [  Ff.  |  degenerate  into  diagonal  ones,  see  Kqs.  (12),  thereby 
permitting  simultaneous  solution  of  the  first  N  equations  of  the  stability  matrix,  Fq.  (17).  Hence 


FA(nf  m)  An  +  FB(n,m)BB  t-Fc(n,  m)Cn  -  0 


and 


FB(n.w)  p  Fc(n,m)  p 

n  r, 

FA(n,  tu)  FA(n,  ui) 


(D.3) 


The  coeffidenta  of  the  remaining  2N  equations,  in  Eq.  (17),  are  reorganized  in  the  form: 


£  [GB(n,a>  +  Gc(n,m)  CJ  =  0 
J  (HB(n,m)  Bn  +  Hc(n,m)  CJ  -  0 


n,m 


**» 


*** 


Gq  (n,  m) 


[(l-v)/2] 

n2/32+[(l-v)/2]t2 


8b>|  +  t2K,(n,m) 


1  a2  fi2  +  [(1  -i-/2)t2J 


1 5„_  +tK,(n.m)  -  t3  K,(n,  m) 

n  mt  1  a 


(D.  4) 


(D.  5) 
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Hr  (n,m) 


1- 


(vn/3) 


Ail 


L  n2jS  +[(l-v)/23tJ  12(RA>  12 (R A)  J 


8Bi|  +  Kj  (q,  m) 


n  i  v  i  ,  i4Ks(n,a) 
-  2t2K  (n,m)  +  — — 3 - r- 

12  (R  A) 


(D.5) 


Note  that  here  the  applied  load  i«  lateral  prat  'ure. 


0.2  FORMULATION  OF  INTEGRALS  KJ(a,n)  FOR  RINGS  WITH  SINUSOIDAL  HEIGHT  VARIATION 
The  variation  of  the  height  of  the  rioga  along  the  generator  ia  aaaumed  to  obey  the  law 

Cu  (x)  =  Cj  ain  /3x  +  c03  (D.  6) 

where  cQ2  ia  the  height  of  a  uniform  ring  distribution,  and  c£  ia  the  maximum  height  of  the  non- 
uniform  part. 

Sinusoidal  variation  ia  analysed  since  both  stiffened  and  unstiffened  shells  buckle  under  lateral 
pressure  with  one  half  sine  wave  along  the  generator. 

The  non-uniform  rings,  Eq.  (D.6),  are  compared  to  uniform  ones  of  the  same  weight.  If  the  width 
of  the  rings  is  the  same  in  both  cases  the  equivalent  height  is 

e2  =  coj  +  cj  2/»r  (D.7) 

Tie ’variation  ratio",  y,  is  defined  by 

y  =  i  ,2/cj  (D*8) 

Substitution  of  (D.7)  and  (D.  8)  into  (D.6)  yields 

c2(x)  =  Cj  [y  +  (rr/2)(l— y)  sin  /3x]  (D.9) 

The  geometrical  parameters  of  the  rings  p2  (x) ,  yj(x)  and  qj2(x)  zre  defined  for  rectangular 
cross-sections  in  Appendix  A.  Substitution  of  Eq.  (D.9)  into  those  definitions  yields 
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2 

pj  (x)  =  (1  -  y  )  (d2/a)  aj  (a  +  ain  fix) 

Xj(x)  =  <[(1  -  v2)/2(R/h)]  (d2/a)  at  (c  +  eia/3x)  [1  +  (a  +  tin  fix)] 

q*2(x)  -  12(1— i/2)(dj/a)  ^j[a2  (a  +  ain  j9x)S]/12}  + [(a  +  aln  /3x)/4]  (1  +  ctj  (a  +  ain  fix)]2"^  fD.  10) 

where 


aj  «  (1  -  y)  (rr/2)(c2/h) 

and  (D.  11) 

a  -  (2y/ rr)/(  1  -  y) 

The  integrala  K^fa.m)  ,  Eqe.  (16)>  are  now  computed  with  aid  of  Eqs.  (D.  10)  and  yield 
Kj  (n,  m)  =  (1  -  v2)(d2/a)  at  [a5e>i  +  Ljfu.m)] 

K2(n,m)  =  r[(l-v2)/2(R/h)](d2/a)aj  {[a+o2  a2 +(81/2)35^ +(l  +  2aaj)  Lj  (n,  m)  - 

K}(n>m)  =  12(l-»/2)(d2/a)aj  ^[(ttJ  a2/3)  +  (a  +  a1)(l+2aa1)/4)  5om +l[(l  +  2aa,)  +a2]/4)  Lj  (n,  m) - 
-(ai/8)(l+2«a1)(S1..|l|f2-S1  ,)  -(«2/12)L2(n,m)}  (D.  12) 


where 


"1.1 


1  n  -  m  »  1 
^  0  a  /  1  or  m  /  1 
*+*. 


and 


»-»  .2 


(  4  ,  [*  +  <-»"  1  [  I  _  1  1 

1  n  [  1  -  (a -m)5  1  -  (0+m)2  J 

-J-J 

~(n+m),i  j 


L2(n,m)  - 


3(1  >(-•■)•**!  [_ 


1, 9  -  (n  -  nip”  9  - 


1  |n-m|  =  2 
0  |n  -n|  t  2 


m  1  •'i 

j  >>} 


£ 


t 

I 
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APPENDIX  -  E 
NON-UNIFORM  STRINGERS 

E.1  REDUCTION  OF  THE  ORDER  OF  STABILITY  DETERMINANT 
For  stringers  only 


V-2  (*)  =  XjM  =  =  *?*2  (*)  =  0  (K-D 

and  therefore  (see  Eqs.  16) 

K.(n,m)  =  0  j  =  1,  2,  3  (K.2) 

Matrices  [GA],  { Gg]  and  [Gc]  degenerate  to  diagonal  ones,  see  Eqs.  1?.  Ihe  second  set 
of  Eqs.  (11)  is  solved  simultaneously 


B  «- 


GA(n,m) 

Gfl  (n.  rr.)  B 


Gc(n.m)  £ 
Gg  (n,  m)  n 


(E.3) 


The  new  coefficients  of  the  starility  determinant  of  Eq.  (17),  are 
2  {FA(n,m)  Aa  +  Fc(n,m)  CB]  =  0 
£  [HA(n,m)  +  Hc(n,oi)  CJ  =  0 


(E.  4) 


where 


=  [(l-i/)/2)(nJ/3J  +  *2)2  s  a.  ,  v 

rA(a,  o»)  - - - — oam+nm/3  lj(n,m) 

v'  +  [(W)/2U30 

f(l-^)/2)n/f(»/n202-t?)  „ 


Fr(n,m)  =  H.(m,  n)  = 


t 2  +  (0  —  1 )  '2j  n2  ft 


onm  -n2a fi  l,(n,m) 


Hr(n, 


c  to,  m/  - 


22  2 

o  j_  .  2\  u  1  \  /m  l  £ 

Ip  Ji  T  !  I  (U-W/4JQ  P 

o 


/  t  /At  2  A^  I 

\\f  ifl 


12(R/b) 


tz  +  [(l-„)/2]n*/3  12(R/h) 


J,.+ 


(E.  5) 


m 


3? 


,  ♦  . 


I3  (n,  m)  +  ts  n  m /3  I^ln.ra) 
12(R/h)? 


The  applied  load  ia  here  axial  compression. 


(E.S) 


E.2  FORMULATION  OF  INTEGRALS  Lin, a)  FOR  STRINGERS  WITH  SINUSOIDAL  HEIGHT 
VARIATION. 

The  assumptions  and  the  details  of  the  analysis  are  similar  to  those  presented  in  the  second 
part  of  Appendix  D  for  rings  with  sinnsoidal  height  variation.  The  subscript  *2"  (for  rings)  is  replaced 
by  subscript  "1"  (for  stringers),  e.g. 

c2  -  c, 

C02  "*  C01 

and  the  appropriate  expressions  are  obtained. 

The  stringer  distribution  corresponding  to  the  ring  distribution  of  Eq.  (D.  9)  is  then 

Cj  (x)  =  c,  [y  +  (l-y)(ir/2)  sin  /3x]  (E.6) 

where  Cj  is  the  equivalent  height  of  stringers  sad  y  is  the  variation  ratio. 

If  one  introduces  the  notations 


a  *  2  y/ n  ( 1  —  y) 

and 

a,  {cj/hi  [<r(l— y)/2] 

the  geometrical  parameters  of  the  stringers  .  defined  in  Appendix  A,  become 

pj(x)  =  (l-t'2)(dl/b)a,  (a  +  sin  /3x) 

y*(x)  =  e (( 1  —  i-'2 )/2 ( R/h ) ]  (d  j/b)  Oj  (c  +  sin  fix)  [1  +at(a  +  »i»  /9*)] 

ijJjfx)  -  ‘.2(1  -  i-~)  (d  j/b)  a,  l(af/12)(a  +sin  /3x)  +  (1/4)  (a  +  sio  /3x) (1  +iij  (a  t-sin  (Sx)]  ( 


(E.7) 


(E.  8) 
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r)*j (x)  -  2(l-v)(djA)  (h/b)  aj  (a+sin/ix) 


(E.8) 


Substitution  of  Eq.  (E.8)  into  Eqs.  (16)  yield*  the  appropriate  integral* 


Ij(n.m)  =  (l-v2)(d1/b)«1M,-+L1(n,m)l 

Ij(n.m)  =  t [(l-i/J)/2(R/h)](dj/b) a,  l(a+«Jfi,  +(a1/2)]S>^  +  (l+2aaJ)Lj (n»m) -(a/4) 

Ij  (n,  m)  =  12(lV)(d1/b)a1{[(l/3)a8a?+(l/4)(a+al)(l  +  2«a1)]6,111  + 

+  Kl/4)  [(1  +  2a Ql)2  +  a2ll  Lj  (n, «)  -  (1/8)  ai  (U 2aai) 2  -S,  ,M«2/12)  L2(n,.4 

I4(n,m)  =  2(l-„)(d,/h)Vb)  a,  (aS,.  +  L, (a, m)]  (E-9) 


Lj  (n,  m)  = 


L2  (a,  n)  = 


L}(n,m)  = 


[!  +  (-!)  3  1 _ 1 _  , 

n  1  -  (n  -  m)2  1  _  (n  +  m)^ 

m  X^| 

3[1  +  (-1)  ]  1  1 

n  ^  9  -  (n  -  m)*  9  -  (*  +  m)2 

+  1 _  1  , 

n  1  -  (a  -  »)*  T-  (.+5» 


(E.  10) 


E.3  FORMULATION  OF  INTEGRALS  l.(n,m)  FOR  STRINGERS  WITH  LINEAR  HEIGHT  VARIATION 

Flf.  i. 


Cj(x)  =  C0X  +  c*(2/5./<r)x  0<x<w/2fi 

c,(x)  =  c01  +c*[2-(Wt)i]  n/2fi<x<*/0 

tf/fi  ■  L/R  (E.  11) 


The  equivalent  height  of  the  stringer*  i* 


c,  =  C01  +  (cj/2) 


I 


(E.  12) 
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Substitution  of  Eq.  (E.  12)  into  Eq.  (E.  11)  yield* 
c,(x)  =  F,  [y  +  2(l-y)(3f3/*)*]  0<x<n/2fi 

Cl(x)  =  Cj  ly  +  2(1— y)  [2-(2/3/ff)x]i  (n/20)  <  x  <  (m/0)  (E.  13) 

where  y  is  the  variation  ratio,  defined  here  by  y=c01/c1 

The  geometrical  parameters  of  the  stringers!  for  0  <  x  <  n/ 2/3  are  obtained  from  Eq.  (E.8)  by 
(a)  permutation  of  sin  /3x  to  sin  (2#/ff)  x  (b)  use  of  now  definitions  for  a  and  Cj ,  see  Eq.  (E.  7) 

a  =  y/2(l— y) 

aj  =  (c  j/h)  [2(1  — y)l  (E.  14) 

The  explicit  expressions  for  integrals  (16)  are  then 
Ij  (n,  m)  =  (l-v2)(dj/b)  aj  [a6>M  +  Lj(n,m)] 

I2(n.m)  =  <((l-i/2)/2(R/h)](dj/b)a,[(a+a2a1)5BIi  +  (1  +  2aaj)  Cj(n,m)  + 

+  Oj  L2(n,  m)] 

(j  (n, m)  =  12(l-^)(djA)  aj  {a[(a2aJ/3)  +  (l  +  2aaj)/4]  Smm  + 

+  ((1  +  2aaj)/2]  C1(n,m)  +  l(«j  (1  +  2acj)]/2l  L2(n»m)  (a?/3)  Cs(n,m)} 

I^n.rn)  =  2(l-v)(d,A)3(b/b)  oJaS,. +L4(n,m)]  (E.  15) 


n  £  in 
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APPENDIX  -  F 

COMBINATION  OF  UNIFORM  STRINGERS  WITH  RINGS  OF  NON  UNIFORM  HEIGHT 

For  uniform  stringers,  the  matrices  of  the  integrals  l^U.m),  Eqs.  (16),  sad  those  whose  ele¬ 
ments  ore  FA(n,m),  Fg(a,m)  sad  Fc(a.m),  nee  Eqs.  (12),  degenerate  into  diagonal  matrices. 
Solution  of  the  first  N  equations  of  (17)  and  simplification  of  the  remaining  matrix  yields 


'  Gb 

Gc 

*  « 

B 

• 

* 

► 

.  Ha 

Hc  . 

C 

L  »  J 

y  =  o 


(F.  1) 


where 


2.2 


Gg  (n, m)  =  ) -  S  +  t2 K,  (n, m) 

(l  +  jijln^/S  +((l-v)/2)t2 


Gj.  (n,  m)  =  Hg  (n,  m)  — 


[(l+»/)/2]a2/S  ) 

- —j - 


tS>iB  +  tlCj  (n,  m)  -  tsK2(n,») 


[  (l+(.1).*p  +[(l-t/)/2]t2  j 

!.  «>  J  ^  .4  m  2  «  2. 2 

,.)  Moi.'g  _  .  »  3  £  > - 

12(R/h)  (1  +  m,)o2/3  +[(l-v)/2]t2 


(t2  +  n2/3  /2) 

~  12(R/M* 


t*  K3  (u,  m) 


SBB  +  Kj(n,m)  -  2t2K  (n, m)  + - « 

12(R/h) 


(F.  2) 


The  applied  load  in  these  equations  is  hydrostatic  pressure. 

For  rings  with  sinusoidal  height  variation  the  values  of  integrals  Kj(n,  m)  are  similar  to  those 
given  in  the  second  section  of  Appendix  D. 


42 


CONCLUSIONS 

The  general  instability  of  cylindrical  ahells  with  non*uniform  stiffeners  was  analysed  for  lateral 
and  hydrostatic  pressure  and  for  axial  compression*  Height  and  width  variations  of  rings  and  stringers 
were  considered  and  considerable  improvements  in  structural  efficiency  were  obtained*  When  most  of 
the  stiffening  is  concentrated  at  the  mid  length  of  the  shell*  an  edge  type  buckling  mode  may  appear 
and  reduce  the  structural  efficiency  drastically*  Hence  it  was  found  advisable  in  many  cases  to  sup* 
plement  the  varying  part  of  the  stiffener  by  a  uniform  part.  The  "variation  ratio"  y  characterizes  the 
extent  of  the  uniform  supplement* 

For  lateral  pressure  loading,  gains  of  70  —80%  in  general  instability  pressure  are  obtained  with 
rings  of  non-uniform  height  at  optimal  values  of  y  .  For  prescribed  loading  conditions  the  gains  re¬ 
sult  in  10  —  25%  weight  saving.  Width  variation  is  less  efficient  and  rings  of  non-uniform  width  yield 
only  15—20%  gains  in  load  but  the  corresponding  weight  savings  remain  at  10  —  20%. 

For  axial  compression,  stringers  v/ith  sinusoidal  and  linear  height  variation  were  studied. 
Sinusoidal  variation  is  found  to  be  superior  to  linear  variation*  and  gains  in  buckling  load  of  30%,  with 
corresponding  weight  savings  of  10% ,  are  obtained.  A  variation  ratio  y  «  0*5  appears  to  be  roughly 
optimal. 

For  hydrostatic  pressure  loading,  a  combination  of  rings  and  stringers  with  a  weight  ratio  10  to 
1,  is  found  to  be  roughly  optimal.  Combinations  of  non-uniform  rings  sod  stringers  under  axial  compres¬ 
sion  and  hydrostatic  pressure  should  be  further  investigated. 

Comparison  of  solution  by  a  straightforward  Galerkin  method  acd  by  the  method  of  "correcting 
coefficients*  indicates  that  in  the  relatively  simple  case  investigated  here,  the  loss  in  accuracy  in  the 
second  method  outweighs  the  numerical  difficulties  of  the  first. 
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(c)  RINGS  WITH  WDTH  VARIATION 

FIG.  2  TYPES  OF  STIFFENER  AREA  VARIATIONS 


FIG.  3  SHELL  WITH  CENTRALLY  CONCENTRATED 
RING 


VARIATION  OF  STIFFENERS  FOR  DIFFERENT  VARIATION  RATIOS 


buckling  shape  with  one  basic  wave  in 


DEFLECTION  SHAPES  FOR  LATERALLY  LOADED  CYLINDRICAL  SHELLS 
STIFFENED  WITH  INSIDE  RINGS  OF  NON-UNIFORM  HEIGHT 
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G.  8  DEFLECTION  SHAPES  FOR  LATERALLY  LOADED  CYLINDRICAL  SHELLS 
STIFFENED  WITH  OUTSIDE  RINGS  OF  NON-UNIFORM  HEIGHT 


buckling  load  far  an  n=1  mode  of  instability 


FIG.  9  TYPICAL  CURVE  OF  PRESSURE  PARAMETER  VERSUS  NUMBER  OF 

CIRCUMFERENTIAL  WAVES 
(LATERAL  PRESSURE) 
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LATERAL  PRESSURE  AND  WEIGHT  SA/ING  FOR  SINUSOBAL  HEIGHT  VARIATION  AND  OPTIMUM  K 


INFLUENCE  OF  SINUSOIDAL  HEIGHT  VARIATION  OF  RMGS  ON  LATERAL  PRESSURE 
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FIG.  15  EFFECT  OF  ECCENTRICITY  ON  THE  STRUCTURAL  EFFICIENCY  OF  RINGS  WITH  SINUSOIDAL 

HEIGHT  VARIATION 
(LATERAL  PRESSURE) 


LATERAL  PRESSURE  AND  WEIGHT  SA/ING  FOR  SINUSOIDAL  WIDTH  VARIATION 


EFFECT  OF  ECCENTRICTTY  ON  THE  RELATIVE  STRUCTURAL  EFFICIENCY  OF 
STRINGERS  WITH  SINUSOIDAL  HEIGHT  VARIATION 
(AXIAL  COMPRESSION) 
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STRUCTURAL  EFRCENCY  OF  TWO  TYPES 
OF  STRINGER  HEIGHT  VARIATION 
(AXIAL  COMPRESSION) 


COMPARISON  OF  STRUCTURAL  EFFICIENCY  OF  LINEAR  AND  SINUSOIDAL  HEIGHT 

VARIATION  OF  STRINGERS 
(AXIAL  COMPRESSION) 


FIG.  23  NFLUENCE  OF  SINUSODAL  HEIGHT  VARIATION  OF  RINGS  ON  HYDROSTATIC 

PRESSURE 


INFLUENCE  OF  SINUSOIDAL  HEIGHT  VARIATION  OF  NSIDE  RINGS  WITH  UNIFORM 

STRINGERS  ON  HYDROSTATIC  PRESSURE 


INFLUENCE  OF  SINUSOIDAL  HEIGHT  VARIATION  OF  OUTSIDE  RINGS  WITH  UNIFORM 
OUTSIDE  STRINGERS  ON  HYDROSTATIC  PRESSURE 


